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We analyse the self-diffusiophoresis of a spherical particle animated by a nonuniform chemical 
reaction at its boundary. We consider two models of solute absorption, one with a specified dis¬ 
tribution of interfacial solute flux, and one where this flux is governed by first-order kinetics with 
a specified distribution of rate constant. We employ a macroscale model where the short-range 
interaction of the solute with the particle boundary is represented by an effective slip condition. 
The solute transport is governed by an advection-diffusion equation. 

We focus upon the singular limit of large Peclet numbers, Pe 1. In the fixed-flux model, the 
excess-solute concentration is confined to a narrow boundary layer. The scaling pertinent to that 
limit allows to decouple the problem governing the solute concentration from the flow field. The 
resulting nonlinear boundary-layer problem is handled using a transformation to stream-function 
coordinates and a subsequent application of Fourier transforms, and is thereby reduced to a non¬ 
linear integral equation governing the interfacial concentration. Its solution provides the requisite 
approximation for the particle velocity, which scales as Pe^ 1 ^ 2 . 

In the fixed-rate model, large Peclet numbers may be realized in different limit processes. We 
consider the case of large swimmers or strong reaction, where the Damkohler number Da is large 
as well, scaling as Pe. In that double limit, where no boundary layer is formed, we obtain a closed- 
form approximation for the particle velocity, expressed as a nonlinear functional of the rate-constant 
distribution; this velocity scales as Pe~ 2 . Both the fixed-flux and fixed-rate asymptotic predictions 
agree with the numerical values provided by computational solutions of the nonlinear transport 
problem. 


I. INTRODUCTION 

When a nonuniform chemical reaction takes place on 
the boundary of a colloidal particle, the interaction of 
the reactants with that boundary results in a flow field 
and a consequent motion of the freely suspended particle. 
When the solute is electrically neutral, its short-range in¬ 
teraction with the particle boundary may be represented 
by an effective slip condition, relating the fluid velocity at 
the outer edge of the interaction layer to the tangential 
gradient of solute concentration [3]. A simple “contin¬ 
uum” model of such a slip-based self-diffusiophoresis was 
provided by Golestanian et al. [7], who for simplicity de¬ 
scribed the chemical reaction by a prescribed distribution 
of solute flux. A different “colloidal” description for this 
problem was also proposed by Ref. [6]; its linkage to the 
continuum approach appears to be controversial Hi- 

A key assumption in the model of Ref. [7] is the ne¬ 
glect of solute advection, resulting in a linear transport 
problem governing the solute concentration. It appears 
that the first systematic analysis of solute advection in 
the continuum description was carried out by Michelin & 
Lauga m , who considered the relatively simple config¬ 
uration of a spherical particle with an axially symmetric 
distribution of chemical reactions. They considered two 
separate models of solute production: in the first, follow- 
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ing Ref. [7], the flux of solute is prescribed along the par¬ 
ticle boundary; in the second, following Ref. U], the rate 
constant associated with a first-order chemical reaction is 
specified there. Starting from the exact microscale model 
together with either a ‘fixed-flux’ or ‘fixed-rate’ boundary 
conditions, and using coarse-graining techniques familiar 
from electrokinetic analyses 18:, 21] , Ref. [13] derived a 
macroscale model where the solute-boundary interaction 
is systematically transformed into an effective slip condi¬ 
tion (cf. Ref. [3]). 

Michelin & Lauga m solved their macroscale problem 
for finite values of the Peclet number Pe using compu¬ 
tations. For the fixed-flux variant of that problem, these 
computations indicate that at large Pe the particle speed 
scales as Pe -1 ^ 3 . While that scaling was already pre¬ 
dicted by [5] , the associated boundary-layer problem has 
not been solved. In this paper we address this singular 
problem using a boundary-layer analysis. This is sup¬ 
plemented by a large-Peclet-number analysis of the com¬ 
parable fixed-rate problem. In both problems we derive 
asymptotic approximations for the particle velocity, and 
compare them with the values obtained from the compu¬ 
tational solution of the respective macroscale models. 


II. PROBLEM FORMULATION 

We consider first the fixed-flux variant of the general 
problem described by Ref. [13] . A sphere of radius a is 
suspended in a solution of an otherwise uniform concen¬ 
tration, say ^oq. Due to surface reaction, the particle 
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emits solute at a prescribed rate srf , which in general 
varies along the particle boundary; it is assumed axially 
symmetric. The interaction energy between the solute 
and the boundary, characterized by the lengthscale A, is 
comparable with the thermal energy kT. Our goal is the 
calculation of the steady-state velocity acquired by the 
particle due to a nonuniform boundary distribution of 
g/. The characteristic value of the excess-solute concen¬ 
tration, relative to the ambient value, is ^ = ||,c/||a/-D, 
where \\g/\\ is some representative value of g/ and D is 
the solute diffusivity. The interaction of the solute with 
the particle boundary is represented as a body force on 
the fluid, whose magnitude is of order kTT?/a; balancing 
this force with the viscous stress in the 0(A)-wide inter¬ 
action layer provides the pertinent velocity scale, namely 

IKII kTX 2 ^ 
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wherein % is a unit vector in the 9 = 0 direction and IA 
is the velocity of the particle relative to the fluid; and 
the condition that the particle is force-free. The flow is 
animated by the slip condition, 

u = MVjC at r = 1, (7) 

wherein M is the diffusio-osmotic slip coefficient and V s 
the surface-gradient operator. For an isotropic interac¬ 
tion potential M is uniform, but may be either positive 
or negative; the scale was chosen by Ref. m to make 
M of unity magnitude, hence 

M = ±1. (8) 

The goal is the calculation of the swimming velocity U 
as a function of Pe and the distribution k(n). Because of 
the standard structure of a slip-driven Stokes-flow prob¬ 
lem, U is obtained using the reciprocal theorem [3 [2t|j, 


r) being the solution viscosity. 

We employ the macroscale formulation of Ref. ’13], ap¬ 
propriate to the limit where the interaction thickness A 
is small compared with a. In that formulation the per¬ 
tinent fields are the flow and excess-solute concentration 
outside the interaction layer. In the dimensionless no¬ 
tation of Ref. m distances are normalized by a, the 
excess-solute concentration by and the velocity scale 
is chosen comparable to ([T]) (see below); the governing 
equations are written in a reference frame moving with 
the particle, where the problem is steady, using spheri¬ 
cal coordinates in which the radial distance r is measured 
from the sphere centre and the zenith angle 9 is measured 
from the symmetry axis: see figure [lja). For convenience 
we actually use the alternative coordinate fj, = cos 9. 

The excess concentration c is governed by the 
advection-diffusion equation, 


Pe u • Vc = V 2 c, 


where 



( 2 ) 

(3) 


is the Peclet number and u the velocity field in the fluid; 
the imposed-flux condition, 


dc 

dr 


k(n) at r = 1, 


(4) 


where k = —g //||,£/|| is an 0(1) scaled activity; and the 
attenuation condition, 


U = • (f u dA (9) 

4tt / r=1 

where dA is a dimensionless areal element (normalized 
by a 2 ). Use of the axisymmetric slip |t| followed by 
integration by parts thus provides the formula 

U = —M J /j,c(r = 1, /u) d/x. (10) 

While the particle speed depends only upon the surface 
concentration of c, one cannot avoid in general the so¬ 
lution of the flow problem: the concentration field is af¬ 
fected by the flow through the advection term in ([2|. 
Since this coupling between the flow and solute concen¬ 
tration is nonlinear, so is the dependence of IA upon both 
Pe and k(fi). 

Note that, in contrast to classical forced-convection 
problems [2], the Peclet number here does not involve an 
externally imposed velocity scale. Rather, it represents 
the intensity of interfacial chemical activity. Specifically, 
substitution of |l]) into (|3j) implies that Pe is of order 

\\^\\kTX 2 a 

r]D 2 [ ’ 

Michclin & Lauga m used computations to analyse self¬ 
propulsion at arbitrary Peclet numbers and regular per¬ 
turbations to obtain analytic approximation for small 
Peclet numbers. In what follows, we employ singular 
perturbations to address the opposite asymptotic limit, 
Pe > 1. 


c —>• 0 as r —► oo. 


(5) III. LARGE-PECLET-NUMBER LIMIT 


The velocity field u is governed by the Stokes equa¬ 
tions; the far-held approach to a uniform velocity of mag¬ 
nitude U (see figure [lja)), 

( 6 ) 


The advection-diffusion equation § suggests that in 
the large-Peclet-number limit the solute concentration is 
dominated by advection, so at leading order 

u ■ Vc = 0. 


u —^ —ilA as r —> oo, 


( 12 ) 
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Since the flow is steady, this implies that c is a constant 
along the streamlines of the leading-order flow. Antici¬ 
pating open streamlines, which originate and end at in¬ 
finity, condition ([5]) then implies 

c = 0. (13) 


The nil result (13) is clearly incompatible with Q. 
This non-uniformity is associated with the singular na¬ 
ture of the large-Peclet-number limit, where the high¬ 
est derivative is effectively multiplied by a small parame¬ 
ter. It suggests that a diffusive boundary layer is formed 
about the boundary r = 1; in that thin layer, advection 
and diffusion of solute are comparable. Note that the 
Stokes flow is coupled to the solute concentration c only 
through its value at r = 1, and is accordingly ‘unaware’ 
of the steep gradient of c near that surface; consequently, 
no boundary layer exists in the flow field. 


A. Boundary-layer formulation 

We denote the thickness of the layer by (>(^C 1). Its 
scaling with Pe is readily obtained using dominant bal¬ 
ances mm- from condition (El) we see that c = 0(<5) in 
the layer. The slip condition Id) then implies the same 
scaling for the tangential velocity-component v , while the 
continuity equation in conjunction with the impermeabil¬ 
ity condition implies that the radial velocity-component 
u is of order S 2 . The advection-diffusion equation 0 
thus yields 5 = Pe 1 ^ 3 . It follows that in the boundary 
layer c is 0(Pe -1 / 3 ); the velocity field u is of the same 
magnitude in the entire fluid domain, and then, given 
([9]), so is also U. 

The above scaling suggests the expansions u = 
Pe^ l ^U + ■ ■ ■ and v = Pe~ 1 ^ 3 V + ■ ■ • for the veloc¬ 
ity components, where the rescaled components are 0(1) 
functions of r and fi. Within the boundary layer we use 
the expansion 


c = Pe~ 1 ^ 3 C + • • • , (14) 

where the rescaled concentration C is an 0(1) function 
of the stretched radial coordinate 


about r = 1, which, in view of ([T6|a), read 

+ V=V\ r=1 + ---. (18) 

Use of the leading-order continuity equation, 

% < 19 > 

allows to express both approximations in terms of V\ r ,. 
Substitution into the advection-diffusion equation (pi) 
and making use of ©) thus yields the parabolic equa¬ 
tion, 




-r|;[(l-^)/Wlg, (20) 

governing the boundary-layer solute concentration. It is 
supplemented by the boundary condition (cf. @) 

3C 

— = k(n) at Y = 0 (21) 

and the matching condition (see ©) 

0 as Y-> 00 . (22) 


The boundary-value problem (20)—(22) is uncoupled to 
that governing the flow, and may be solved indepen¬ 
dently. In view of (10) there is actually no need to solve 
for the flow: once C is known, the particle velocity is 
obtained from (10) and (14). Note however that with 
f(fj) being related to C through ©> equation ( p0| ) is a 
nonlinear one. 

The preceding boundary-layer problem appears simi¬ 
lar to that governing the nutrient concentration about a 
squirming sphere mum- In the feeding problem, how¬ 
ever, the slip distribution is prescribed, and is accordingly 
independent of the Peclet number; this results in a linear 
problem with a boundary-layer thickness that scales as 
the —1/2 power of that number. 


B. Boundary-layer analysis 


Y = Pe 1/3 {r- 1) (15) 


and /r. In terms of the rescaled components, condition 
0 becomes 

{7 = 0, V = —M( 1 - /x 2 ) 1/2 /({t) at r = 1, 

(16a, b) 

where 


f(p) 


def dC 
dji 


y=o 


(17) 


Following the analysis of a similar boundary-layer 
problem HZ!. our approach in confronting (p0[)-(p2| is to 
calculate C as if were a prescribed quantity. The so¬ 
lution of the linear system (20)—(22) in conjunction with 
© then provides a nonlinear equation governing C at 
Y = 0. Given (10) and (14), this interfacial distribution 
is all that we need to obtain the swimming velocity. 

Our scheme is appropriate for the class of problems 
where the boundary layer does not detach; given ([7| this 
necessitates that on the boundary C is strictly increasing 
or decreasing function of /r; with no loss of generality we 
assume the latter, 


Since U and V vary on an 0(1) scale they are approx¬ 
imated in the boundary layer using Taylor expansions 


f(p) < 0 . 


(23) 
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Given © and the vanishing of C at the outer edge of 
the boundary layer, it appears plausible that for (23) to 
be satisfied k)fi) must be strictly increasing. 


The parabolic equation ( [20] ) does not appear amenable 
to integral transform. We accordingly shift to stream- 
function coordinates [TOj, replacing the independent vari¬ 
ables by (^, 7 ?), with 


* = —Y{1 - /x 2 )/(/x) 


(24) 


(which is proportional to the Stokes stream-function near 
r = 1.) Substituting this change of variables into (20) we 
find that (7(>^,/x) satisfies the equation 




while condition (21) becomes 

k{ii) 


dC 


(i-^ 2 )/(m) 


at if - = 0; 


(25) 


(26) 


given (23), the decay condition (22) now applies as 'S 
oo: 


C —► 0 as if - —> oo. 


(27) 


The boundary-value problem (25)—(27) is naturally 


solved using a Fourier cosine transform, defined as 


CW) = 


,\ 1/2 

0 . 


C{^, n) cos(w!F) dw. (28) 


Application of this transform to (25) (27) results in the 
first-order equation 


^ + Mu 2 (l - /x 2 )/(/i)C = M 0 


1/2 


k{n). (29) 


The complementary solution of (291 is A(w)e 
where 


x(aO = fiq)il-q 2 )dq; 
J o 


(30) 


given ( |23|, x(/i) is strictly decreasing. A particular inte¬ 
gral ofl29| is 


M 


7 fme- M “ 2 ^-xWd£. (31) 


The value of /io (and then A) is determined by the 
condition that C —> 0 as w —> oo regardless of the value 
of ii. We now show that, upon choosing /io = A/, this 
condition yields A = 0. Indeed, since the function x is 
strictly decreasing this choice results in a negative power 
in the exponent appearing in (31), whereby the particular 


solution vanishes as lo —> oo. Given (23), the sign of y(/i) 
is opposite to that of fi, implying that the exponent in the 
complementary solution diverges as w -> oo for a finite 


interval of /i-values (0 < /i < 1 for M = 1; — 1 < /i < 0 
for M = —1). Since C must vanish in that limit for all 
H, it follows that A = 0. 


We conclude that C is provided by (31) with /x 0 = 


M. Applying the inverse transform, interchanging the 
order of integrations, and evaluation at = 0 we obtain, 
upon substitution of © and ( |30| ), the following equation 
governing C on the particle boundary (denoted C(/i) for 
brevity): 


C'(M) = 4 / 

J A 


M 


m d£ 


lMj\l- q ^dq 


(32) 


d q 


Since (7(/i) is strictly decreasing, the local slip is directed 
in the ±9 direction for M = ±1. Given (10), this im¬ 


plies that the particle moves in the ±i direction. In the 
particle-fixed reference frame, the direction of the inci¬ 
dent flow is then In view of the parabolic nature 
of the boundary-layer problem, the value of C at a cer¬ 
tain value of fi is only affected by the distribution of C 
upstream. For M = 1 this upstream distribution corre¬ 
sponds to the interval (/x, 1), while for M = —1 it cor¬ 
responds to (—1,/x). Equation (32) is indeed consistent 


with this information-propagation property. 

Equation ( |32[ constitutes an integral equation gov¬ 
erning Cin) [cfT [T]. Upon expanding C(/x) in Legendre 
polynomials, C(/x) = A n P n ill), it is transformed 

into a nonlinear algebraic system governing the coeffi¬ 
cients {A n }. This system may be solved using an it¬ 
erative scheme when C(/x) is strictly decreasing. Once 
solved, the particle velocity is obtained using (10) and 
©, yielding 


U k, —MAiFe" 1/3 . 
3 


(33) 


For a given distribution fc(/x) one needs in general to solve 
( [32] ) separately for M = ±1, resulting in two different 
distributions — say C±iii). In the case where fc(/x) is an 
odd function, however, it is readily verified from (32) that 


C-iii) = ~C + i~ii). Given (10), the resulting particle 
velocities are opposite in sign. This, of course, is obvious 
by symmetry. 


The asymptotic prediction (33) can be compared with 
a direct computational solution of equations (@) ob¬ 
tained using a spectral decomposition in the azimuthal 
direction and a stretched radial grid (see Ref. [T3]| for 
more details). As an example we consider the simple dis¬ 
tribution kifi) = ii. Since this function is odd, it is suffi¬ 
cient to solve (32) for M = 1: we then get from (10) that 

0.57Pe _1 'A 


A\ = —0.86 whereby ( ].‘{.'l[ ) gives U « 0.57Pe _i/,:i . This 
result is in excellent agreement with the full numerical so¬ 
lution: see figure 0b). The solute-concentration maps, 
shown in the insets, illustrate the transition from an es¬ 
sentially fore-aft symmetric distribution at weak convec¬ 
tion to a boundary-layer structure at strong convection. 

Note that the use of the macroscale model of Ref. [T3| 
within a boundary layer of dimensional thickness a5 im- 
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FIG. 1: (a) Schematic, (b) Particle velocity U for the fixed-flux model, with k(p) = p and Mm 1. Dashed line: large-Pe 
approximation (331; symbols: computational results. The insets depict the solute-concentration distribution (scaled by the 
maximum surface concentration) as provided by the numerical solutions for Pe = 1 and Pe = 100. 


plies that the underlying limit of that model, A -C a, 
should be refined to A <C aS. Moreover, the very deriva¬ 
tion of the effective conditions © and ([TJ) actually re¬ 
quires (A /a)Pe < 1 and (A /a) 2 Pe <C 1, respectively [ Kil . 
With the boundary-layer thickness S being Pe -1 / 3 it fol¬ 
lows that the stringent condition which must be satis¬ 
fied is A /a -C Pe^ 1 . Ref. [T3] estimated that A/a is 
between 10~ 5 and 10 -3 , implying that the computations 
performed for the largest Peclet number employed herein 
(Pe = 100) still fall within the validity domain of the 
model. 


IV. FIXED-RATE MODEL 


We now consider the fixed-rate model, where the rate- 
constant of a first-order kinetic relation is specified 
on the particle boundary (and assumed axially symmet¬ 
ric). Since a characteristic rate of solute absorption is 
ll^ll^oo, where ||J^|| is a characteristic norm of J^, the 
velocity scale ([!]) is replaced with 

imi^oo^A 2 

ijD 1 1 

The Peclet number is now of order (cf. ©) 
a\\Jf\\^kT\ 2 

r^D2 ' 

The limit of short-range interaction is again governed 
by© ©, except that © is replaced by 

dc 

= fc(/i)(l + Dac) at r = 1 


(35) 


dr 


wherein 


Da = 


ape 1 

D 


(36) 


(37) 


is the Damkohler number and now k = Jff /||J^||, which 
is again an 0(1) distribution; since the first-order kinetic 
describes solute adsorbing onto the boundary, k is here 
non-negative. 

When addressing here the limit of large Peclet numbers 
care must be exercised. Given (35), one would typically 


envision large Pe due to fast reaction (large ||JP||) or 
large swimmers (large a). Since however both Pe and 
Da are linear in a||JP||, the proper limit to consider is 
that where Da becomes O(Pe). 

With a large Peclet number, one may naively expect 
the topology of the fixed-flux analysis, namely the trivial 
solution (13) with a non-zero excess concentration con¬ 


fined to a thin boundary layer. A dominant-balance in¬ 
spection reveals however that no boundary-layer struc¬ 
ture is compatible with condition (361. This apparent 


paradox is readily resolved by noting that, if 
u = 0(Pe~ 1 ) or asymptotically smaller, 


(38) 


the advective term does not dominate © Approxima¬ 
tion (12) is then rendered invalid, as is then (13); thus, 


no boundary layer is realized. 

We therefore proceed under the a priori assumption 


(38). The scaling of c and u are accordingly obtained by 


considering the dominant balances of (36) in the absence 


of a boundary layer. As it turns out, the only consistent 
balance is that between the last two terms in (36), im¬ 
plying that c is 0(Da~ 1 ). Note however that this is not 
the velocity scaling, since the present balance implies a 
uniform leading-order value of c (namely —1/Da) on the 
boundary, so the slip is triggered by the leading-order 
correction to the 0(Da~ 1 ) concentration. It follows from 
(36) that this correction is of order Da ~ 2 , thus providing 


the velocity scaling. With Da = 0(Pe ), this is indeed 
compatible with (381. 
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With the fluid velocity being 0(Pe~ 2 ) it becomes 
evident that the leading-order transport is (counter¬ 
intuitively) unaffected by advection. Following the pre¬ 
ceding arguments we postulate the expansions 

c = Da~ 1 Ci+Da~ 2 C 2 +- • • , u = Da~ 2 u 2 +- ■ ■ , U = Da~ 

(39) 


At 0(1), condition (36) yields 


Ci = — 1 at r = 1. 


(40) 


From 0 we find that the leading-order concentration ci 
is harmonic. The solution that satisfies (40) and decays 
at infinity is the monopole 


ci = - 

r 


At 0(Da *), condition (36) reads 


dc\ 

—— = kui)C 2 at r — 1. 
or 


(41) 


(42) 


Substitution of (41) yields c 2 = 1 /fc(/t) at r = 1. Note 
that there is no need to solve for c 2 in the fluid domain, 
where it is governed by an advection-diffusion equation. 
Indeed, using (10) we readily obtain 


Un = -M 


/id/t 

-l Kv)' 


(43) 


In figure [2] we compare the asymptotic approximation 
U w Da~ 2 U 2 with the computational solution. This is 
done for M = 1 and the rate-constant distribution = 
1 + /z/2, for which (43) yields U 2 = 2 In 9 — 4 « 0.3944. 


The computations were performed for Da = Pe. At large 
Da they indeed agree with the asymptotic approxima¬ 
tion. The solute-concentration maps, shown in the in¬ 
sets, illustrate the transition at large Da to the radially 
symmetric distribution (41) 


V. CONCLUDING REMARKS 

We have analysed self propulsion of a chemically reac¬ 
tive particle for large Peclet numbers. In the fixed-flux 
model, the excess-solute concentration is localised in a 
narrow boundary layer. Use of boundary-layer approxi¬ 
mations reduces the coupled problem governing the non¬ 
linear solute transport and fluid motion to the solution 
of an integral equation governing the interfacial solute 
concentration. In the fixed-rate model, we focused upon 
the limit of large swimmers or strong reaction, where the 
Damkohler number becomes comparable to the Peclet 
number. In that problem no boundary layer occurs, and 
solute advection actually diminishes with increasing Pe. 

The respective scaling of the dimensionless velocity 


with Pe 


-1/3 


and Da z implies that the velocities 0 
and (34), characteristic of the flow at moderate Pe , are 



FIG. 2: Particle velocity U for the fixed-rate model, with 
fc(/i) = 1 + /x/2 and M = 1. Line: large-Pe approximation; 
symbols: computational results, performed for Da = Pe. The 
insets depict the solute-concentration distribution (scaled by 
the maximum surface concentration) as provided by the nu¬ 
merical solutions for Da = 1 and Da = 50. Note the approach 
to a spherically symmetric distribution with increasing Pe 
(cf. Ell). 


of ([T]) and |TT] ) in the fixed-flux case and (34) and (371 in 
the fixed-rate case, we find instead the respective velocity 
scales (cf. Ref. 01) 


||^|| 2 / 3 (UT) 2 / 3 A 4 / 3 c tf 00 kT\ 2 D 


^2/3£)l/3 a l/3 


?ya 2 ||jr|| 


(44) 


representing a transition from a size-independent velocity 
to ones that decrease with particle size. With the generic 
velocity scales (!]) and ( [34] ) being non-representative at 
large Pe, it follows that the dimensionless number Pe 
itself, as provided by © and ( |35| ), does not constitute 
a genuine Peclet number, but simply a measure of the 
surface activity. This is reminiscent of electrokinetic phe¬ 
nomena at strong applied fields mTnm or due to im¬ 
posed flows 03], where the generic Peclet number defined 
by Ref. (II] no longer represents the relative magnitudes 
of advection and diffusion. 

The present work suggests two future directions. The 
first is a different analysis of the fixed-reaction-rate 
model, appropriate to the case where Pe becomes large 
due to large solute molecules (small D). In that scenario, 
where Da is 0(Pe 1 ' 2 ), the boundary-layer thickness is 
Pe ~ 3 / 2 with an identical scaling for the velocity field (so 
(38) is satisfied). The bound ary-lay er problem is then 
similar to that formulated in (IIIA except that (21) is 


no longer representative at large Pe. Rather, making use 


replaced by an inhomogeneous Robin condition. 

The other direction involves the case where the chem¬ 
ical reactions produce ions, rather than neutral species. 
In that case the interaction layer is the Debye diffuse- 
charge layer. This problem is fundamentally different, 
as the interaction potential (namely the electric poten¬ 
tial) is itself coupled to the ionic concentrations through 
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Poisson’s equation. This “auto-electrophoresis” problem 
therefore falls into the realm of electrokinetics [II] . The 
macroscale description of that problem, appropriate to 
the limit of thin double layers, was developed by Ref. [21] ■ 
It differs from the macroscale description of Ref. m in 
several fundamental aspects, the important one being the 
need to solve for the nonlinearly-coupled transport of two 
fields (ionic concentration and electric potential). Solute 
advection then constitutes only one of several nonlinear¬ 
ities inherent in the problem. The macroscale model of 
Ref. [22j was solved in that paper using a linearization 
scheme, appropriate to the case of a nearly homogeneous 


particle. It is desirable to extend this solution with nu¬ 
merical computations, similar to those of Ref. [13], as 
well as a large-Peclet-number asymptotic analysis, com¬ 
parable to that appearing in the present paper. 
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